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The interaction between arbitrary dielectric heterostructures is studied within the framework of 
a recently developed dielectric contrast perturbation theory. It is shown that periodically patterned 
dielectric or metallic structures lead to oscillatory lateral Casimir-Lifshitz forces, as well as modu- 
lations in the normal force as they are displaced with respect to one another. The strength of these 
oscillatory contributions increases with decreasing gap size and increasing contrast in the dielectric 
properties of the materials used in the heterostructures. 



OS 

o 
o 

(N 



PACS numbers: 05.40.-a, 81.07.-b, 03.70.+k, 77.22.-d 



3 



> 

so 



I. INTRODUCTION 



In light of the ongoing miniaturization of mechanical devices and the recent developments in Casimir-Lifshitz 
interactions [H, [HSU HI, there has been some recent interest in the effect of these interactions between the components 
of small mechanical devices @. Since these interaction are particularly strong at small distances, it will be interesting 
to know how they can be utilized for designing novel mechanical systems that could work without physical contact 
\ and could potentially help solve the wear problem Q ■ 

In the past few years there have been a surge of interest in developing techniques that can be used to study the 
Oh' Casimir-Lifshitz interaction in non-ideal geometries, including geometry perturbation theories [9|. [i p| , semiclassical 
+L> . [HI and classical ray optics fl2| approximations, multiple scattering and multipole expansions jl3l. Il4 [ill llH \vi\, 
world-line method [l8| and exact numerical diagonalization methods [iH , [20| , as well as the numerical Green's function 
calculation method [2l|. These methods have been used in studying the Casimir force in a variety of different 
geometries, which have improved significantly our understanding of the nontrivial geometry dependence of this effect. 

The effect of non-ideal geometry has been shown to lead to a number interesting effects. For example, it has 
been suggested that corrugated surfaces opposite one another can experience an oscillatory lateral Casimir force 
Q, which was subsequently observed experimentally Q. A recent experiment probing the normal Casimir force 
between a smooth surface and surface with tall rectangular corrugations also revealed further evidence on the non- 
additive nature of the Casimir force Here, we study the Casimir-Lifshitz interaction between arbitrary dielectric 
(f-) , heterostructures within the framework of a recently developed formalism [2^, [27| • We derive a closed form expression 
£SJ ' for the Casimir-Lifshitz energy between two dielectric heterostructures (such as the example depicted in Fig. [TJ up 
to the second order in the perturbation theory and show that a coherent coupling between the different modes of the 
spectrum of the dielectric pattern takes place across the gap. As a special example, we consider unidirectional periodic 
^\ ' heterostructures (see Fig. [I} and calculate the lateral and normal Casimir-Lifshitz force between them within the 
same order in the perturbation theory. We find that coupling between modes with identical wavevectors of the pattern 
structures between the different objects can lead to modulations in the normal force and can give rise to oscillatory 
• i-h , later forces, reminiscent of the lateral Casimir force that appears due to coupling between geometrical features such 
^ 1 as corrugations [3, Q . 

This paper is organized as follows. Section |TT] sketches the dielectric contrast perturbation theory, and Sec. IIIII 
elaborates on how it can be used for dielectric heterostructures giving closed form expressions for the second order 



term in the perturbation theory. Section IIVI gives the results for the lateral and normal Casimir-Lifshitz force for a 
number of choices of materials, and Sec. [V] contains some discussions and concluding remarks. 



II. THEORETICAL FORMULATION 



To calculate the Casimir-Lifshitz interaction we need to quantize the electromagnetic field in a background that 
includes the dielectric or metallic objects that modify the quantum fluctuations of the field. Describing a general 
assortment of dielectric and metallic objects in space via a frequency dependent dielectric profile e(w, r), we can write 
a general expression for the Casimir-Lifshitz energy as [13] 

^ = ^|°°§trln[/C y (C;r,r')], (1) 
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FIG. 1: Schematic representation of two identical semi-infinite and periodic objects made of intercalated layers of high and 
low dielectric functions, occupying the fractions of / and 1 — /, respectively. Here H is the separation between them, a is a 
dimensionless lateral displacement, and A is the wavelength of the periodic structure. 



where 



ICij — 



c 2 

-^e(iC,r)S tj + didj - d k d k Si 



<5 3 (r-r') 



(2) 



We can consider the dielectric function as e(i(, r) = 1 + Se(i(,r), and expand Eq. CD ) in powers of the dielectric 
contrast. A similar approach has been the subject of a few recent studies [23. l23l uM. l25l. l2o| . 

The expansion leads to the decomposition of ICij into a diagonal part /Co,r/ > corresponding to the empty space, and 
a perturbation part SfCij , namely 



1Cij (C; q, q') = JC ,ij (C, q)(27r) 3 5 3 (q + q') + <5%(C; q, q'), 



(3) 



where 



c 2 



(4) 



and 



c 2 

<5£ij(C;q,q') = — <Me«,q + q' 



(5) 



This yields an expansion 



where 



°° C— I 1 )™ -1 

trln[/C] = trln[/C ] + V - - tr[(/C7 1 <$£)"], 



^«(C,q) 



(6) 



(7) 



The first term is the vacuum energy in the absence of the objects, and the terms in the series take account of their 
effect in a perturbative scheme. The n-th order term in Eq. ([6|) takes on the explicit form 



tr[(JCo ^JC)"] = J 



(27T)3 (2tt) 3 [g+g(l)2]...[g + g(„)2] 



L«5e(iC,-q«+q (2) ) ■ • • 5e(i{, -q<"> +q«), 

(8) 
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which involves the Fourier transform of the dielectric contrast profile. Going to real space, we can rewrite the energy 
of the system as [27j 



E = h 
where 



dC A (-l) r 



zi \ " 

2tt ^ 



n=l 



d 3 ri ■ ■ ■ d 3 r n A ili2 (ri - r 2 ) • • • A inil (r„ - n 



5e«,ri) 



l + i<ye(iC,ri) 



5e(iC,T n ) 
l + |<5e(iC,r„) 



A,(C,r) 



2 p -Cr/c 



c 2 Anr 



Sa 1 



c c 
1 + 3— + 3 



£r ' ~ Q 2 r 2 



(9) 



(10) 



We now use this formulation to study the Casimir-Lifshitz interaction between structures with inhomogeneous or 
patterned dielectric properties. 



III. DIELECTRIC HETEROSTRUCTURES 



Let us now consider a configuration similar to the one depicted in Fig. [T] namely two dielectric heterostructures 
that are placed parallel to each other at a separation H. Using the definition r = (x, z), the dielectric profile can be 
written as 



' £u«,x), 

e(*C,r) = " 1, 

, £d«,x), 



H 



< z < +oo, 



< z < f 1 



-oo < z < 



—H 



using the labels u and d for the "up" and "down" bodies respectively. 

To keep the calculations tractable, we now focus on the second order term in the series expansion in Eq. 
such two semi-infinite bodies, the second order interaction term between the bodies can be written as 



E 2 = -- 



2n 2 c 2 



d( C 2 / <fxd 2 x' 



d 2 Q 

(27T)* 



Q-(x-x') 



£(Q) 



1 + i5e„«,x) 



Se d (iC,x') 



l + Ue d (i£,x>) 



for any lateral dielectric function profile, where 



£{Q) 



dp 



[2p 4 - 2p 2 + 1] 



[4p 2 + (cQ/0- 



i3/2 



(11) 



For 



(12) 



(13) 



and 6e U)d (iC,x) = e„ ld «,x) - 1. 

We now focus on the specific example of unidirectional periodic structures as depicted in Fig. [TJ which is made of 
subsequent layers of materials with relatively high and low dielectric functions. We can use the periodic properties of 
the dielectrics and write them in Fourier series expansion. As Fig. [1] shows, we can define the dielectric profile of the 
d-object as 



r e i (io , 



sA < x < + sX, 



e d (iC,x) = < 



e h «) , + sX < x < Q + sX, 



(14) 



[ ei (»C) , ^ + sA < x < 



sX, 



where s is an integer number. We define the Fourier series as 

Se d (iC,x) 



1 + kSe d (i(,x) 



i2irmx / X 



where 



C m {iO 



sin mnf 



m— — oo 



Se h (*C) 



Set (*0 



1 + \5e h (»C) 1 + i<5ez (iC) 



(15) 



(16) 
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for m/0, and 

Co^O = f 



Se h (*0 



1 + Ue h (*C) 



+ (1-/) 



1 + |tfei (i<) 



(17) 



We can find the corresponding expansion for the u-object by changing x — > x + aA. 

Using the Fourier series expansion, one can find the Casimir-Lifshitz energy between two dielectric heterostructures 
as depicted in Fig. [1] [up to second order in the Clausius-Mossotti expansion of Eq. fl5J)] as 

E pp=-^I1 /"^^(^r) C(*C)cos(2™), (18) 

where the prime on the summation sign indicates that the m = term is counted with half the weight, and the pp 
index means the energy calculated for the plate-plate geometry. This result shows that similar to the case of two 
corrugated surfaces, two patterned dielectric heterostructures also couple to each other at the leading order when 
the two wavelengths of the modulations are equal Moreover, higher harmonics contribute to the Casimir-Lifshitz 
energy with exponentially decaying contributions, such that at large separations only the fundamental mode (lowest 
harmonic) will survive [19[ . 

IV. THE NORMAL AND LATERAL FORCES 

We now use Eq. (JT5J) to calculate the normal and lateral forces between different types of dielectric and metallic 
heterostructures. We look at three different types of materials as examples, namely, gold, silicon, and air/vacuum, 
and consider layered materials made of gold-silicon, silicon-air, and gold-air. We describe the dielectric function of 

gold using a plasma model, namely, €{iQ = 1 + -^j, where ui p is the plasma frequency, which is given as w p (Au) = 



w2 



1.37 x 10 16 rad/s [28|. For silicon we use the Drude-Lorentz form e{iQ = 1-1 — - — ^—5, where w p (Si) = 3.3 wo(Si) and 

— C + w o 

wo (Si) = 6.6 x 10 15 rad/s [28(. Finally, for air/vacuum we use e(iQ = 1. 

Due to difficulties in keeping the surfaces of the objects parallel to each other, most experiments are performed 
in plate-sphere geometry. To perform the calculation of the forces for the plate-sphere configuration, we can use the 
Derjaguin Approximation (29j . where we replace one of the semi-infinite objects with a planar surface with a sphere 
with radius R. The approximation is valid provided that the radius of sphere is much larger than the distance between 
the dielectric heterostructures, namely, R 3> H. Using this approximation we can find the normal force between a 
semi-infinite dielectric heterostructure and a sphere of the same material composition as [29] 

F«T = 2nR (^f) . (19) 



, A , 

Using this result, we can find the Casimir-Lifshitz energy for plate-sphere configuration as 



IH 

which we can now use to calculate the lateral Casimir force as 

1 dR, 
A da 



E PS = dH' F™ r (H>), (20) 



^? = -T^- (21) 



Substituting Eqs. US} and |20]) into Eq. ([21]). it reads 



The above equations are the basis of the results that will be presented below. 

Figures [2^,-c show the normal Casimir-Lifshitz force between two unidirectional (layered) dielectric heterostructures 
as shown in Fig. [T]when laterally displaced with respect to one another by aA. It corresponds to the symmetric case 
with / = 0.5, and the corrugation wavelength of A — 1 /im. Three different compositions of gold-silicon, silicon-air, 
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FIG. 2: Normal Casimir-Lifshitz force between layered dielectric heterostructures as shown in Fig. \T\ in the plate-sphere 
geometry, corresponding to (a) gold-silicon, (b) silicon-air, and (c) gold-air, with / = 0.5, and to (d) gold-silicon, (e) silicon-air, 
and (f) gold-air, with / = 0.2. The numerical value of corrugation wavelength used is A = 1 /im. Different curves correspond 
to different gap sizes of H — 100 nm, H = 300 nm, and H = 600 nm. The forces are normalized to Fp S that corresponds to 
the normal force when the laterally-averaged dielectric profile is used. 



and gold-air are considered each at three different gap sizes of H = 100 nm, H = 300 nm, and H = 600 nm. The 
normal forces are normalized using the normal force F® s that corresponds to the Casimir-Lifshitz force calculated 
within the same scheme but with laterally averaged dielectric profile, which corresponds to the m = term in the 
expansion in Eq. (|18p . The normal force is found to oscillate as a function of the lateral displacement, having the 
maximum value when the regions of high dielectric constant from both sides are exactly opposite one another, and 
the minimum value when in the staggered configuration where regions of higher dielectric constant face regions of 
lower dielectric constant. The amplitude of the oscillations increases by decreasing the gap size, and the effect is 
progressively stronger when the contrast between the dielectric properties of the two regions is more pronounced, 
with a maximum relative change of 0.7 % for gold-silicon, 7 % for silicon-air, and 65 % for gold-air, at the closest 
separation of H = 100 nm. 

In Figs. [2ji-f the normal Casimir-Lifshitz forces between the same types of structures as above are presented, for the 
asymmetric case of / = 0.2. One can see two noticeable differences with the symmetric case. First, the oscillations are 
now asymmetric, as enforced by the asymmetry of the dielectric profile, although the asymmetry weakens as the gaps 
size increases and eventually disappears — i.e. the oscillations become symmetric and harmonic — at sufficiently large 
separations. This is consistent with the picture that different harmonics of the dielectric contrast profile in Eq. (fl8)) 
couple with each other via an exponential terms that decays with the corresponding wavelengths of each harmonic 
and as a result any asymmetry caused by higher harmonics will die out at large gap sizes. The second new feature 
is the significant enhancement of the amplitude of the oscillatory behavior as a function of the lateral displacement. 
While it is still the case that this amplitude increases with increasing contrast between the dielectric properties of the 
two materials used in the layered structure, the maximum relative change is 0.4 % for gold-silicon, 6 % for silicon-air, 
and 200 % for gold-air, at the closest separation of H — 100 nm. 

The lateral Casimir-Lifshitz forces for the same layered structures as above are shown in Figs. [3^-c for the symmetric 
case with / = 0.5. In this case, we have assumed R = 180 fim and A = 1 fim. Similar to the previous study, three 
different compositions of gold-silicon, silicon-air, and gold-air are considered each at three different gap sizes of H = 100 
nm, H = 200 nm, and H = 400 nm. The lateral force is found to oscillate as a function of the lateral displacement, 
reminiscent of the lateral Casimir force that is induced by geometrical corrugations 0, Q . The shape of the oscillatory 
function approaches a sinusoidal behavior as the gap size increases, consistent with the fact that higher harmonics do 
not contribute to the force in that limit as also seen in geometrical lateral Casimir effect [l!|. The amplitude of the 
oscillations increases by decreasing the gap size as well as the contrast between the dielectric properties of the two 
regions. Numerically, we find an amplitude of 0.5 pN for gold-silicon, 8 pN for silicon-air, and 12 pN for gold-air, at 
the closest separation of H = 100 nm. 

Figure[3Ji-f show the lateral Casimir-Lifshitz forces between the same types of structures as above, for the asymmetric 
case of / = 0.2. Similarly, the profiles of the lateral force are noticeably asymmetric, with the asymmetry weakening 
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FIG. 3: Lateral Casimir-Lifshitz force between layered dielectric heterostructures as shown in Fig. Q] in the plate-sphere 
geometry, corresponding to (a) gold-silicon, (b) silicon-air, and (c) gold-air, with / = 0.5, and to (d) gold-silicon, (e) silicon-air, 
and (f) gold-air, with / = 0.2. The numerical values used in these graphs are A = 1 /im and R — 180 /^m. Different curves 
correspond to different gap sizes of H = 100 nm, H = 200 nm, and H — 400 nm. 



as the gap size is increased and the shape of the profile approaches that of a sinusoidal function (single harmonic). 
We also see comparatively more significant enhancement of the amplitude of the oscillatory behavior as a function of 
the lateral displacement. The amplitude of the oscillations is found as 0.3 pN for gold-silicon, 5 pN for silicon-air, 
and 7 pN for gold-air, at the closest separation of H = 100 nm. 



V. DISCUSSION 



In this paper, we have proposed a mechanism by which it is possible to create a lateral Casimir-Lifshitz force as 
well as controlled modulations in the normal Casimir-Lifshitz force without geometrical corrugations. A coupling 
similar to what exists in the case of corrugated surfaces gives rise to these oscillatory forces, namely identical modes 
of the dielectric patterns couple across the gap to generate a macroscopic coherence in the fluctuations. The generic 
features of these oscillatory forces are very similar to those of the forces caused by corrugations; the effect is stronger 
and involves more harmonics at closer separations, while it weakens and only involves the lowest mode of the pattern 
in the dielectric contrast at larger separations. 

While the difference in the dielectric properties of the materials controls the general strength of the above results, 
comparison between Fig. [2] and [3] shows that the modulations in the normal force are more strongly affected by the 
contrast in the dielectric properties. The choice of air/vacuum as one component also allows us to make predictions 
about geometrical features with large corrugation amplitudes, which provides an approximation scheme for the non- 
perturbative geometrical regime. 

In the present calculations we have only used the second order terms in the dielectric contrast perturbative series. 
Higher order terms shown in Eq. ||5J) will introduce coupling between different modes of the dielectric pattern in a 
systematic way, as imposed by the overall conservation of the sum of all wavevectors (momenta) . While the present 
is aimed at showing in terms of tractable calculations, one can in principle carry out the calculation of the Casimir- 
Lifshitz interaction in such dielectric heterostructures using numerical diagonalization methods [201 ] . 

Controlled interactions between dielectric heterostructures with smooth outer surfaces could be very useful in 
practical applications because it will help avoid the complications of bringing surfaces with geometrical protrusions 
close to each other while avoiding contact between them and controlling their separations. Moreover, it is much easier 
to pattern dielectric properties of materials in a controlled way than it is to shape them with the high precision that 
is needed for Casimir effect type experiments. 



7 



Acknowledgments 

The authors thank the ESF Research Network CASIMIR for providing excellent opportunities for discussion on the 
Casimir effect and related topics. This work was supported by EPSRC under Grant EP/E024076/1. 



[1] H.B.G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 

[2] E.M. Lifshitz, Sov. Phys. JETP 2, 73 (1956); I.E. Dzyaloshinskii, E.M. Lifshitz, and L.P Pitaevskii, Adv. Phys. 10, 165 
(1961). 

[3] S.K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); U. Mohideen and A. Roy Phys. Rev. Lett. 81, 4549 (1998); A. Roy and 
U. Mohideen, Phys. Rev. Lett. 82, 4380 (1999); H.B. Chan et al, Science 291, 1941 (2001); G. Bressi et al, Phys. Rev. 
Lett. 88, 041804 (2002); R.S. Decca et al Phys. Rev. Lett. 91, 050402 (2003); D.E. Krause, Phys. Rev. Lett. 98^ 050403 
(2007). 

[4] F. Chen, U. Mohideen, GL. Klimchitskaya, and V.M. Mostepanenko, Phys. Rev. Lett. 88, 101801 (2002); Phys. Rev. A 

66, 032113 (2002); 
[5] H.B. Chan et al, Phys. Rev. Lett. 101, 030401 (2008). 

[6] F.M. Serry et al, J. Microelectromech. Syst. 4, 193 (1995); E. Buks and M.L. Roukes, Phys. Rev. B 63, 033402 (2001); 
H.B. Chan et al, Phys. Rev. Lett. 87, 211801 (2001); G. Palasantzas and J.Th.M. De Hosson, Phys. Rev. B 72, 121409 
(R) (2005). 

[7] A. Ashourvan, M.F. Miri, and R. Golestanian, Phys. Rev. Lett. 98, 140801 (2007); Phys. Rev. E 75, 040103 (R) (2007); 
T. Emig, Phys. Rev. Lett. 98, 160801 (2007); M. Miri and R. Golestanian, Appl. Phys. Lett. 92, 113103 (2008); F.C. 
Lombardo et al, J. Phys. A 41 164009 (2008); I. Cavero-Pelaez et al, Phys. Rev. D 78 065019 (2008). 
[8] R. Golestanian and M. Kardar, Phys. Rev. Lett. 78, 3421 (1997); Phys. Rev. A 58, 1713 (1998); M. Kardar and R. 

Golestanian, Rev. Mod. Phys. 71, 1233 (1999). 
[9] T. Emig, A. Hanke, R. Golestanian, and M. Kardar, Phys. Rev. Lett. 87, 260402 (2001); Phys. Rev. A 67, 022114 (2003). 
[10] P.A. Maia Neto, A. Lambrecht, and S. Reynaud, Phys. Rev. A 72, 012115 (2005); R.B. Rodrigues, P.A. Maia Neto, A. 

Lambrecht, and S. Reynaud, Phys. Rev. Lett. 96, 100402 (2006). 
[11] M. Schaden and L. Spruch, Phys. Rev. A 58, 935 (1998). 
[12] R.L. Jaffe and A. Scardicchio, Phys. Rev. Lett. 92, 070402 (2004). 

[13] R. Balian and B. Duplantier, Ann. Phys. (New York) 104, 300 (1977); 112, 165 (1978). 
[14] O. Kenneth and I. Klich, Phys. Rev. Lett. 97, 160401 (2006) 
[15] R. Golestanian, Phys. Rev. E 62, 5242 (2000). 

[16] T. Emig, N. Graham, R.L. Jaffe, and M. Kardar, Phys. Rev. Lett. 99, 170403 (2007). 

[17] T. Emig, N. Graham, R.L. Jaffe, and M. Kardar, Phys. Rev. A 79, 054901 (2009). 

[18] H. Gies and K. Klingmuller, Phys. Rev. Lett. 96, 220401 (2006); Phys. Rev. Lett. 97, 220405 (2006). 

[19] T. Emig, Europhys. Lett. 62, 466 (2003); R. Biischer and T. Emig, Phys. Rev. Lett. 94, 133901 (2005). 

[20] A. Lambrecht and V.N. Marachevsky, Phys. Rev. Lett. 101, 160403 (2008). 

[21] A. Rodriguez, M. Ibanescu, D. Iannuzzi, F. Capasso, J.D. Joannopoulos, and S.G Johnson, Phys. Rev. Lett. 99, 080401 
(2007); A. Rodriguez, M. Ibanescu, D. Iannuzzi, J.D. Joannopoulos, and S.G. Johnson, Phys. Rev. A 76, 032106 (2007); 
A. Rodriguez, J.D. Joannopoulos, and S.G. Johnson, Phys. Rev. A 77, 062107 (2008). 

[22] G. Barton, J. Phys. A 34, 4083 (2001). 

[23] R. Golestanian, Phys. Rev. Lett. 95, 230601 (2005). 

[24] S.Y. Buhmann and D.-G. Welsch, Appl. Phys. B 82, 2, 189 (2006). 

[25] G. Veble and R. Podgornik, Eur. Phys. J. E 23, 275279 (2007). 

[26] K.A. Milton, P. Parashar, and J. Wagner, Phy s. Rev. Lett. 101 , 160402 (2008). 

[27] R. Golestanian, Phys. Rev. A, in press (2009) |arXiv:0905.T046] . 

[28] E. D. Palik, Handbook of Optical Constants of Solids, edited by E. D. Palik (Academic Press, New York, 1985). 
[29] J.N. Israelachvili, Intermolecular and Surface Forces (Academic, London, 1992). 



Casimir-Lifshitz Interaction between Dielectric Heterostructures 



Arash Azari, Himadri S. Samanta, and Ramin Golestanian 
Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, United Kingdom 

(Dated: July 14, 2009) 

The interaction between arbitrary dielectric heterostructures is studied within the framework of 
a recently developed dielectric contrast perturbation theory. It is shown that periodically patterned 
dielectric or metallic structures lead to oscillatory lateral Casimir-Lifshitz forces, as well as modu- 
lations in the normal force as they are displaced with respect to one another. The strength of these 
oscillatory contributions increases with decreasing gap size and increasing contrast in the dielectric 
properties of the materials used in the heterostructures. 



OS 

o 
o 

(N 



PACS numbers: 05.40.-a, 81.07.-b, 03.70.+k, 77.22.-d 



3 



> 

so 



I. INTRODUCTION 



In light of the ongoing miniaturization of mechanical devices and the recent developments in Casimir-Lifshitz 
interactions [?????], there has been some recent interest in the effect of these interactions between the 
components of small mechanical devices [? ]. Since these interaction are particularly strong at small distances, it 
will be interesting to know how they can be utilized for designing novel mechanical systems that could work without 
| \ physical contact and could potentially help solve the wear problem [? ] . 

In the past few years there have been a surge of interest in developing techniques that can be used to study the 
Oh' Casimir-Lifshitz interaction in non-ideal geometries, including geometry perturbation theories [???], semiclassical 
+i» . [? ] and classical ray optics [? ] approximations, multiple scattering and multipole expansions [? ? ? ? ? ], world- line 
method [? ] and exact numerical diagonalization methods [? ? ] , as well as the numerical Green's function calculation 
method [? ]. These methods have been used in studying the Casimir force in a variety of different geometries, which 
have improved significantly our understanding of the nontrivial geometry dependence of this effect. 

The effect of non-ideal geometry has been shown to lead to a number interesting effects. For example, it has 
been suggested that corrugated surfaces opposite one another can experience an oscillatory lateral Casimir force [? 
], which was subsequently observed experimentally [? ]. A recent experiment probing the normal Casimir force 
between a smooth surface and surface with tall rectangular corrugations also revealed further evidence on the non- 
additive nature of the Casimir force [? ]. Here, we study the Casimir-Lifshitz interaction between arbitrary dielectric 
, heterostructures within the framework of a recently developed formalism [? ? ] . We derive a closed form expression 
for the Casimir-Lifshitz energy between two dielectric heterostructures (such as the example depicted in Fig. ??) up 

^ | \ to the second order in the perturbation theory and show that a coherent coupling between the different modes of 

the spectrum of the dielectric pattern takes place across the gap. As a special example, we consider unidirectional 
Q\ ' periodic heterostructures (see Fig. ??) and calculate the lateral and normal Casimir-Lifshitz force between them 
within the same order in the perturbation theory. We find that coupling between modes with identical wavevectors 
of the pattern structures between the different objects can lead to modulations in the normal force and can give rise 
• i— i , to oscillatory later forces, reminiscent of the lateral Casimir force that appears due to coupling between geometrical 
^ ■ features such as corrugations [? ? ]. 

This paper is organized as follows. Section ?? sketches the dielectric contrast perturbation theory, and Sec. ?? 
elaborates on how it can be used for dielectric heterostructures giving closed form expressions for the second order 
term in the perturbation theory. Section ?? gives the results for the lateral and normal Casimir-Lifshitz force for a 
number of choices of materials, and Sec. ?? contains some discussions and concluding remarks. 



II. THEORETICAL FORMULATION 



To calculate the Casimir-Lifshitz interaction we need to quantize the electromagnetic field in a background that 
includes the dielectric or metallic objects that modify the quantum fluctuations of the field. Describing a general 
assortment of dielectric and metallic objects in space via a frequency dependent dielectric profile e(w, r), we can write 
a general expression for the Casimir-Lifshitz energy as [? ] 

^ = ft^°°§trln[/C y (C;r,r')], (1) 
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FIG. 1: Schematic representation of two identical semi-infinite and periodic objects made of intercalated layers of high and 
low dielectric functions, occupying the fractions of / and 1 — /, respectively. Here H is the separation between them, a is a 
dimensionless lateral displacement, and A is the wavelength of the periodic structure. 



where 



c 2 

-^e(i(,r)5ij + d l d j - d k d k 5ij 



S 3 (r-r r ). 



(2) 



We can consider the dielectric function as e(i£,r) = 1 + Se(i(, r), and expand Eq. (??) in powers of the dielectric 
contrast. A similar approach has been the subject of a few recent studies [?????]. 

The expansion leads to the decomposition of /Qj into a diagonal part )Co,ij, corresponding to the empty space, and 
a perturbation part SK-ij , namely 



where 



fCij (C; q, q') - tC ,ij (C, q)(27r) 3 5 3 (q + q') + 5JC l3 (C; q, q'), 



c 2 



(3) 



(4) 



and 



c 2 

(5/Cy(C;q,q') = -^S ll Si(iC 1 q + q')- 



(5) 



This yields an expansion 



trln[/C] =trIn[/C ] + 5^ 

n=l 



(-1)' 



tr[(/C^J/C) n ] 



(6) 



where 



^«(C,q) = 



£1 



(7) 



The first term is the vacuum energy in the absence of the objects, and the terms in the series take account of their 
effect in a pcrturbativc scheme. The n-th order term in Eq. (??) takes on the explicit form 



tr[(/C -^/C)"] = 



(2tt) 3 (2^) 3 



rC 2 



+ 9 (1)2 ] • • • [£ + <Z ( " )2 ] 



^C,-q (1) +q (2) ) • • -^C-q^+qW), 

(8) 
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which involves the Fourier transform of the dielectric contrast profile. Going to real space, we can rewrite the energy 
of the system as [? ] 



roc ^> oo (_1)n-l r 

E = h -± V i — / d 3 n ■ ■ ■ d 3 r n A llt2 (n - r 2 ) • • • A inh (r„ - n) 



Se(i(, ri 



l + i(5e(iC,r„) 



where 



A,(C,r) 



2 p-Cr/c 



C 2 e 



c 2 A-rrr 



1 + 



£r ' £ 2 r 2 



c c 
1 + 3— + 3 



(9) 



(10) 



We now use this formulation to study the Casimir-Lifshitz interaction between structures with inhomogeneous or 
patterned dielectric properties. 



III. DIELECTRIC HETEROSTRUCTURES 



Let us now consider a configuration similar to the one depicted in Fig. ??, namely two dielectric heterostructures 
that are placed parallel to each other at a separation H. Using the definition r = (x, z), the dielectric profile can be 
written as 



e(«C,r) 



1, 



H 



< z < +oo, 



2 ^ Z ^ 2 ' 



(11) 



-oo < z < 



-H 



using the labels u and d for the "up" and "down" bodies respectively. 

To keep the calculations tractable, we now focus on the second order term in the series expansion in Eq. (??). For 
such two semi-infinite bodies, the second order interaction term between the bodies can be written as 



Ed 



2tt 2 c 2 ./„ 



d( C 2 / d 2 xd 2 x' 



d 2 Q 

(2n) 2 



iQ-(x— x') 



1 + Ue u (iC,x) 



5e d (i(,x' 



1 + Ue d (i(,x>) 



for any lateral dielectric function profile, where 

£{Q)-- 



dp [2p 4 -2p 2 + l]^ e _^V4p 2 +(cQ/C) 2 . 



[V + (cQ/0 



213/2 



(12) 



(13) 



and 5e u , d {i(,x.) = e M , d «,x) - 1. 

We now focus on the specific example of unidirectional periodic structures as depicted in Fig. ??, which is made 
of subsequent layers of materials with relatively high and low dielectric functions. We can use the periodic properties 
of the dielectrics and write them in Fourier series expansion. As Fig. ?? shows, we can define the dielectric profile of 
the d-object as 



' ei(i(), -*+ s \<x<-Q + s\, 

e d «, x) = I e h «) , + sX < x < Q + sA, 

k ei «) , Q + sX < x < | + sA, 
where s is an integer number. We define the Fourier series as 

5e d (i(,x) 



(14) 



1 + \5e d (iC,x) 



where 



C m (i() = 



sin rrmf 



= C m {iQ^ mx '\ 



Sei (*C) 



5e h (iQ 
1 + Ue h (*C) 1 + \5ei (*C) 



(15) 



(16) 



